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On  Shortest  Paths  in  Polyhedral  Spaces 

Micha  Sharin*'  and  Amir  Schorr 

School  of  Mathonatkal  Sciences 
Tel  Aviv  Umversity 

ABSTRACT 

We  consider  the  problem  of  computing  the  shortest  path  between  two  points 
in  two-  or  three-dimensional  sp3ace  bounded  bv  polyhedral  surfaces.  In  the  2-D 
case  the  problem  is  easily  solved  in  time  C?(nnog  n).  In  the  general  3-D  case 
the  problem  is  quite  hard  to  solve,  and  is  not  even  discrete;  we  present  a  doubly- 
exponential  procedure  for  solving  the  discrete  subproblem  of  determining  the 
sequence  of  boundary  edges  through  which  the  shortest  path  passes.  Finally  we 
consider  a  favorable  special  case  of  the  3-D  shortest  path  problem,  namely  that  of 
finding  the  shortest  path  between  two  points  along  the  surface  of  a  convex 
polyhedron,  and  solve  it  in  time  0(n    log  n). 

1.   Introductioo 

The  problem  of  tlnding  the  shortest  path  between  two  points  in  Euclidean  space  bounded  by 
a  finite  collection  of  polyhedral  obstacles  is  a  special  case  of  the  more  general  problem  of  planning 
optimal  collision-free  paths  for  a  given  robot  system  (here  we  treat  the  robot  as  a  single  moving 
point). 

In  two  dimensional  space  the  problem  is  easy  to  solve,  because  the  shortest  pjath  between  two 
given  points  must  be  a  polygonal  line  whose  vertices  are  comers  of  the  given  polygonal  obstacles, 
so  that  the  problem  can  be  immediately  reduced  to  a  discrete  graph  searching,  and  can  be  solved 
in  time  0(n  log  n),  where  n  is  the  number  of  obstacle  comers.  This  two-dimensional  problem 
has  been  considered  by  Lozano-Perez  and  Wesley  [LW],  and  later  also  by  Lee  and  Prcparata  [LP]. 
In  some  special  cases,  considerably  more  efficient  algorithms  exist.  For  example,  if  the  free  space 
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within  which  the  shortest  path  is  sought  is  the  interior  of  a  simple  pxjlygon,  then  the  shortest  path 
can  be  found  in  time  0{n  log  n)  [LP],  [Ch].  As  another  example,  if  all  the  barriers  are  straight 
segments  parallel  to  each  other,  then  the  problem  can  again  be  solved  in  time  0{n  log  n)  [LP]  (a 
similar  favorable  case,  but  with  a  somewhat  different  solution,  is  discussed  below).  Another 
favorable  case  is  noted  by  Torapa  [To],  where  the  obstacles  are  all  convex  and  aligned  in  a  certain 
manner  along  a  straight  line. 

In  three-dimensional  space  the  problem  becomes  much  harder.  In  this  case  the  shortest  path 
between  two  given  points  can  also  be  shown  to  be  a  polygonal  line,  but  all  we  can  say  about  its 
vertices  is  that  they  lie  on  edges  of  the  given  polyhedral  obstacles.  Thus  the  problem  is  by  no 
means  discrete.  Even  if  one  knew  the  sequence  of  obstacle  edges  through  which  the  desired 
shortest  path  passes,  the  calculation  of  the  points  of  contact  of  the  path  with  these  edges  requires 
solution  of  high-degree  algebraic  equations,  which  must  be  accomplished  either  by  numerical 
approximate  methods,  or  by  precise,  but  very  inefficient,  symbolic  algebraic  calculations.  Even 
the  calculation  of  the  sequence  of  obstacle  edges  through  which  the  shortest  path  passes  seems  to 
be  very  difficult,  and  we  do  not  know  of  better  than  doubly-exponential-time  algorithms  for  this 
subproblem.  Papadimitriou  [Pa]  has  recently  presented  an  approximating  algorithm  for  the 
general  three-dimensional  polyhedral  shortest  path  problem,  which  runs  in  pseudo  polynomial 
time. 

However,  in  certain  special  cases  the  3-D  problem  is  not  so  hard  to  solve.  We  will  consider  the 
case  of  finding  the  shortest  path  between  two  points  along  the  surface  of  a  convex  polyhedron,  a 
problem  that  has  been  suggested  originally  by  KE.  Dudcney  as  a  mathematical  puzzle  in  1903 
(see  [Ga,  p.  36];  in  his  original  formulation  a  spider  has  to  aawl  along  the  surface  of  a  cube  to 
reach  a  fly  in  the  shortest  possible  manner).  We  present  an  0{n^  log  n)  algorithm  for  this 
problem,  which  exploits  the  special  structure  of  "geodesic"  paths  along  the  surface  of  a  convex 
polyhedron.  Our  technique  has  later  been  extended  by  CRourke,  Suri  and  Booth  [OSB]  to  find 
shortest  paths  along  the  surface  of  a  nonconvcx  polyhedron,  and  also  been  improved  by  Mount 
[Mo]. 
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The  paper  is  organized  as  follows.  Section  2  presents  a  straightforward  solution  to  the  2-D 
problem;  Section  3  discusses  the  general  3-D  problem,  develops  the  relevant  theory,  and  presents 
a  doubly-exponential  algorithm  for  solving  the  problem.  Finally,  Section  4  analyzes  the  problem  of 
shortest  paths  along  a  convex  polyhedron  and  an  algorithm  for  solving  this  problem  is  presented  in 
Section  5.  ... 

2.  The  Two-dimeasioaal  Case 

If  one  assumes  that  the  solid  obstacles  are  all  orthogonal  prisms  whose  heights  are  all 
parallel  to,  say,  the  z-axis,  and  are  enclosed  between  a  floor  and  a  ceiling,  and  if  the  two  points 
X,  Y,  between  which  an  optimal  path  is  sought,  are  assumed  both  to  have  the  same  height,  then 
the  3—D  case  of  our  problem  may  be  reduced  to  the  2— D  case,  since  then  the  optimal  path  will 
lie  entirely  in  a  horizontal  plane  containing  X  and  Y,  and  the  given  solid  obstacles  will  intersect 
this  plane  in  a  collection  of  polygonal  obstacles.   Let  us  therefore  assume  that  V  is  a  dosed  two- 
dimensional  region  bounded  by  a  collection  of  polygonal  walls  and  other  polygonal  obstacles,  and 
let  X,  Y  be  two  points  in  V.  The  2-D  shortest  path  problem  is  then  to  find  a  (Eudidean)  shortest 
path  between  X  and  Y  which  is  wholly  contained  in  V.  This  problem,  and  several  spedal  cases  of 
it,  has  already  been  considered  by  several  other  people  ([LW],  [LP],  [Ch],  [To]).  The  general 
approach  used  below  has  been  suggested  by  Lozano-Perez  and  Wesley  [LW],  but  without  any 
explidt  analysis  of  its  complexity.  Most  of  the  other  relevant  work  has  involved  spedal  cases  of 
the  problem,  where  more  effident  solutions  than  the  one  presented  below  exist.   To  simplify  the 
following  discussion,  we  assume  that  each  comer  of  the  boundary  of  V  is  inddent  to  just  two 
edges  (although  the  method  described  below  will  also  apply  in  cases  where  this  assumption  does 
not  hold).   It  is  easily  seen  that  the  shortest  route  fromX  to  Y  whidi  is  wholly  contained  in  V  is  a 
polygonal  path  connecting  X  to  y  whose  intermediate  comers  are  all  vertices  of  the  polygonal 
walls  and  obstades  bounding  V.   Furthermore,  a  vertex  p  will  foUow  another  vertex  q  in  such  a 
route  only  if  p  and  q  are  visible  from  each  other  (in  V),  i.e.  only  if  the  straight  segment  joining 
p  and  q  lies  wholly  in  V.  See  Fig.  2.1  for  illustration  of  the  coix»pts  just  discussed. 


/ 


Fig.  2.1.  The  shortest  path  problem  in  2-D  polygonal  space. 
Thus,  to  solve  the  2— D  version  of  oiir  problem,  we  first  construct  a  visibility  graph  VG, 
whose  nodes  are  X,  Y,  and  the  vertices  of  the  boundary  of  V,  and  each  of  whose  edges  connects  a 
pair  of  vertices  visible  from  each  other,  and  has  length  equal  to  the  distance  between  these 
vertices.  Then  we  search  through  VG  to  find  the  shortest  path  in  VG  from  X  loY.  This  yields  a 
polynomial-time  discrete  algorithm  which  can  be  implemented  to  run  in  time  0{nrlogn).  For  the 
sake  of  completeness,  we  sketch  here  such  a  straightforward  implementation. 

Let  N  denote  the  set  containing  X,Y,  and  all  the  vertices  of  the  boundary  of  V  (W(V)  for 
short).  For  each  u  ^  N,  ■wt  will  find  all  points  in  N  visible  from  u  by  using  the  following 
algorithm  which  is  based  upon  a  'plane-sweeping'  technique  similar  to  those  used  in  [Sh],  [NP]. 
Specifically,  for  each  orientation  6  let  /  =  Iq  denote  a  ray  extending  from  u  at  orientation  9.  Let 
Af(9)  denote  the  list  of  all  edges  in  bd{V)  whose  interior  intersects  /,  sorted  in  increasing 
distance  of  these  intersections  from  u.  We  maintain  M(Q)  as  a  2-3  tree,  and  can  easily  compute 
an  initial  value  of  Af  for  some  starting  orientation  Oq  in  time  0{n  log  n). 
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We  next  rotate  the  ray  /  in  clockwise  direction  about  u,  and  update  the  value  of  Af  (9)  each 
time  /  crosses  an  orientation  6  at  which  M  changes,  either  due  to  addition  of  new  edges  into  M  or 
deletion  of  edges  from  M,  or  both.  It  is  dear  that  such  changes  can  occur  only  at  orientations  6 
of  /  at  which  /  passes  through  another  comer  in  N .  We  thus  iterate  through  these  critical 
orientations  in  clockwise  order  (Lncluding  also  the  orientations  of  the  rays  connecting  u  \.oX,  Y). 
Let  6„  be  such  an  orientation  at  which  /  passes  through  p  €  N .  Then  M  will  change  at  6^  so  that 
the  edges  incident  to  p  and  presently  in  Af  are  deleted  from  M ,  and  the  other  edges  incident  to  p 
are  inserted  into  Af.  However,  if  the  ray  /  passes  simultaneously  through  several  points  p  €  N , 
then  we  have  to  remove  from  Af  all  edges  incident  to  any  of  these  comers,  and  to  add  to  Af  all 
other  such  incident  edges,  except  for  edges  which  lie  on  /,  which  are  ignored  in  this  updating 
process. 

At  each  such  orientation  O-  we  also  check  whether  the  first  edge  in  Af  has  changed.  If  it  did  then 
that  comer  p  on  I  incident  to  this  edge  is  visible  from  u,  and  we  update  VG  by  adding  to  it  the 
edge  {u,p).  Otherwise  none  of  these  p's  is  visible  from  u.  If  p  is  X  or  y,  then  it  will  be  visible 
from  u  if  and  only  if  it  precedes  along  /  the  first  edge  in  Af . 

Iterating  in  this  manner  through  all  critical  orientations  6„,  we  find  all  points  in  N  visible 
from  u  in  time  0{n  log  n).  HeiKe,  if  we  repeat  this  procedure  for  each  M  €  N,  we  can 
construct  the  visibility  graph  VG  in  time  0(n  log  n). 

Having  constructed  the  visibility  graph,  we  can  then  search  through  it  to  find  the  shortest 
path  in  VG  connecting  X  and  Y,  using  Dijkstra's  algorithm  (see  e.g.  [AHU]),  which  will  run  in 
time  O(n^). 

Improving  the  Effldency  of  the  Algorithm. 

We  next  discuss  an  alternative  approach  to  the  2-D  shortest  path  problem  which  leads  in 
some  special  cases  to  more  efficient  algorithms.   We  begin  by  the  following  observation: 

Definition:  For  each  point  Z  €  V,  let  it(Z)  denote  a  shortest  path  from  X  to  Z  tlirough  V. 

Lemma  2.1:  Let  Z^,  Z2  €  V.  Then  either  it(Zi)  and  'It(Z2)  do  not  intersect  each  other  or,  if 


they  do  intersect  at  some  last  point  Z,  then  Z  must  be  a  comer  in  W(V),  and  the  lengths  of  the 
initial  portions  of  both  paths  between  X  and  Z  are  equal. 

Proof:  K  tt(Zj)  and  Tr(Z2)  meet  at  a  point  Z,  then  the  lengths  of  their  initial  portions  up  to  Z 
must  be  equal,  or  else  we  could  replace  the  longer  such  initial  portion  by  the  shorter  one,  and  so 
shorten  the  length  of  one  of  these  paths.  The  same  argument  also  implies  that  Z  must  be  a 
comer,  because  otherwise,  after  replacing  one  initial  poitioii  by  the  other  one,  we  would  obtain  a 
shortest  path  to  one  of  the  points  Zi,Z2  which  is  polygonal  and  sias  a  comer  at  an  interior  point 
of  V,  contradicting  the  basic  properties  of  such  shortest  ^laths  nocei  above.   Q.E.D. 

We  can  use  this  observation  to  obtain  sxi  0{n  log  n)  algorithm  for  finding  the  shortest  path 
in  the  following  special  case  (which  is  similar  to  the  second  special  case  considered  in  [LP]): 
Suppose  that  the  boundary  of  V  consists  of  k  v-in;'tkal  lines,  denoted  l^,  .  .  .  ,lf^,  each  of  which 
contains  several  point  apertures  through  which  one  tan  cr>jss  fiom  one  side  of  the  line  to  the 
other.  The  passage  is  blocked  however  at  all  other  points  on  these  barriers.  Suppose  further  that 
X  lies  to  the  left  of  all  these  barriers  and  that  Y  lies  to  the  right  of  all  of  them.  Note  that  the 
shortest  path  from  X  to  y  through  V  must  pass  through  exactly  one  aperture  at  each  of  the  lines 
ll,  .  .  .  .Ij^.  Hence  it  can  be  found  in  0(n^  time,  using  a  standard  dynamic  programming 
approach.    However,  using  Lemma  2.1,  we  can  improve  this  procedure  as  follows  (a  similar 

divide-and-conquer  approach  has  been  used  by  Reif  [Re]  for  a  different  problem  involving  planar 

k 
networks):  Suppose  that  the  barrier  Ij  has  rij  apertures,  ^=1,  .  .  .  ,Jb  (so  that  ^  ";  ~  ^)-   We 

will  process  the  barriers  from  left  to  right.  For  each  barrier  Ij  we  will  compute  for  each  aperture 
Z  €  Ij  the  length  d{Z)  of  the  shortest  path  from  X  to  Z  through  V,  and  also  the  aperture 
p(Z)  €  Ij-i  through  which  the  shortest  path  passes  just  before  reaching  Z.  For  each  Z  €  /^  we 
put  J(Z)  =  lXZ|andp(Z)  =  X. 

Suppose  that  these  maps  have  already  been  computed  for  all  apertures  lying  on  Ij-i-  Let  Z  be  the 
median  of  all  apertures  along  /,.  We  compute  d(Z)  and  p(Z)  by  trying  to  pass  the  path  tt(Z) 
through  each  of  the  apertures  in  Ij^i  (as  in  the  standard  dynamic  programming  approach).  Let 
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W  =  p(Z),  and  assume  that  W  is  the  m-th  highest  point  along  /y_j.  Since,  by  Lemma  2.1, 
shortest  paths  can  be  assumed  not  to  cross  each  other,  it  follows  that  we  can  partition  the  problem 
into  two  subproblems:  First  find  the  shortest  paths  leading  to  the  highest  half  of  the  apertures 
along  /,-,  using  only  the  highest  m  apertures  along  //-^  as  possible  predecessors  along  such  paths, 
and  then  repeat  this  procedure  for  the  lowest  half  of  the  apertiu-es  on  Ij,  using  this  time  only  the 
lowest  n,_i  — m+ 1  apertures  along  Ij..^  as  possible  predecessors. 

Repeating  this  procedure  recursively,  it  is  easily  seen  that  it  will  find  (correctly)  the  shortest 
paths  to  all  apertures  on  /;  in  time  O  ((«,-.- j  +  n^)  log  rij).  Hence  if  wc  iterate  in  this  manner 
through  all  barriers  Ij,  we  obtain  an  0{n  log  n)  algorithm  for  finding  the  desired  shortest  path 
froraX  to  y. 

The  special  case  just  considered  has  led  to  a  fhvorsbis  algorithm  because  of  the  regxilar 
structure  of  shortest  paths  in  this  case  Other  spedni  cases  have  been  considered  in  \LP],  [Qi], 
[To].  In  the  general  case  shortest  paths  may  behave  less  regtilarly,  although  they  still  do  not 
intersect  each  other.  In  fact,  it  is  e-asily  seen  that,  for  a  given  starting  point  X,  the  set  A 
containing  all  the  corners  of  bd{V)  and  X,  can  be  arranged  in  a  tree  T  with  X  as  the  root,  such 
that  each  comer  u  is  the  son  of  a  point  v  if  the  last  straight  segment  on  'It(u)  is  va.  Moreover, 
for  each  u  ^  A  let  r(u)  denote  the  set  of  all  points  y  €  V  for  which  the  last  segment  on  'rr(l') 
is  uy.  Note  that  r(M)  is  nonempty  only  if  the  angle  within  V  between  the  last  segment  vu  on 
it(m)  and  one  of  the  edges  e  of  bd(y)  incident  to  u  is  greater  than  180  degrees.  In  this  case 
r(u)  is  contained  in  the  wedge  formed  between  e  and  the  straight  ray  continuing  vu  past  u.  It  is 
also  easy  to  show,  by  techniques  similar  to  those  used  to  analyze  Voronoi  diagrama  (cf.  [Sh], 
[FAV])  that  the  boundary  arcs  between  adjacent  regions  r(M)  are  all  straight  or  hyperbolic  arcs, 
and  that  there  are  at  most  0{n)  such  arcs. 

In  summary,  the  collection  of  shortest  paths  within  V  from  some  fixed  starting  point  X  can 
be  characterized  by  a  combinatorial  structure  whose  size  is  0{n),  and  it  therefore  seems  likely 
that  faster  than  quadratic  algorithms  for  its  construction  should  exist.  There  exist  some  other 
special  cases  where  this  is  indeed  the  case.  For  example,  if  V  is  the  interior  of  a  simple  polygon. 
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then  shortest  paths  within  V  can  be  computed  in  time  0(n  log  n)  (cf.  [LP],  [Ch]);  another 
favorable  case  has  been  noted  by  Torapa  [To]  in  connection  with  wire  routing  problems  in  VLSI. 
However,  for  general  polygonal  regions  V  the  problem  of  computing  shortest  paths  within  V  in 
fp:.er  than  (9(n^log  n)  time  is  still  open. 

3.   The  Three  Dimensional  Case 

The  situation  becomes  much  more  complicated  when  we  pass  to  the  3-dimensional  version  of 
the  problem.  Here  it  is  easy  to  check  that  the  shortest  path  from  X  lo  Y  consists  of  a  polygonal 
path  whose  vertices  (except  for  X  and  K)  lie  on  some  of  the  edges  of  W(V).  The  problem 
therefore  is  not  immediately  seen  to  be  discrete,  since  there  seems  to  be  a  oontinuura  of  potential 
paths  to  be  considered.  We  will  see  however  that  the  problem  can  be  discretized,  and  develop  an 
algorithm  for  finding  the  shortest  path  which  runs  in  doubly  exponential  time  in  the  number  of 
wall  edges. 

We  win  find  it  useful  to  regard  the  wall  edges  as  open  segments,  so  that  the  wall  comers  are 
disjoint  from  the  edges.  The  collection  of  wall  edges  and  comers  will  sometimes  be  refered  to  as 
"wall  objects". 

We  will  first  consider  the  following  subproblem:  Given  a  sequence  ^  =  (^i,  .  •  .  ,^„)  of  wall 
objects,  find  the  shortest  path  tt  from  X  to  Y  constrained  to  pass  through  each  of  the  objects 
tl»  •  •  •  »^n  '"  ^s  order,  assuming  that  no  other  constraint  is  being  imposed  on  the  path.  Let 
IT  consist  of  the  segments  ttq,  .  .  .  ,1T„. 

Lemma  3.1:  For  each  /  =  1,  .  .  .  ,n,  if  ^,  is  a  wall  edge  (rather  than  a  wall  comer)  then  the 
angles  that  tt,-]^  and  ir,  subtend  at  ^,  are  equal. 

Proof:  This  is  well  known,  but  to  see  this  take  the  two  planes  formed  by  ir/-i  and  ^,  and  by  it, 
and  ^,,  and  "unfold"  them  about  ^^  so  as  to  make  them  coincident,  with  tt^-i  and  it,  lying  on 
different  sides  of  ^,.  Since  tt  is  the  shortest  possible,  the  two  segments  tTj-i  and  it,  must  be 
coUinear  in  this  common  plane,  and  the  angles  which  they  form  with  ^^  must  therefore  be  equal. 
Q.E.D. 
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Let  us  temporarily  assume  that  all  wall  objects  ^i,  .  .  .  ,^„  axe  segments,  rather  than 
points.  An  initial  attempt  to  exploit  Lemma  3.1  in  finding  the  required  shortest  path  it,  is  to 
determine  the  point  Z^  of  contact  of  tt  with  Ij.  This  point  determines  the  angle  at  which  ttq 
meets  ^i,  which  must  be  equal  to  the  angle  at  which  iTi  leaves  ^j.  Knowing  this  angle,  we  can 
determine  the  point(s)  Z2  at  which  tti  will  intersect  ^21  ^^^  continue  to  proceed  in  this  manner 
through  all  segments  in  4.  A  correct  choice  for  Z^  is  one  in  which  the  final  angle  at  which  ir^ 
must  leave  ^„  from  their  point  of  contact  Z„,  is  eq\ial  to  the  angle  subtended  at  |„  by  Z„Y. 

However,  this  approach  is  problematical  in  the  sense  that  each  time  we  try  to  extend  the 
path  TT  from  a  point  Z,  on  ^,  to  ^,+1,  there  can  exist  two  points  Z,  Z'  on  this  edge  for  which  the 
two  angles  subtended  at  ^,  by  the  segments  Z^Z  and  Z^'  are  both  equal  to  the  angle  subtended 
at  ^i  by  IT,-!.  Thus,  even  if  we  fix  the  first  point  of  contaa  Z^,  the  number  of  paths  that  may 
arise  in  the  way  described  above  may  be  exponential  in  n.  To  overcome  this  difficulty,  we  will 
make  use  of  the  following  observation: 

Lemma  3.2:  Let  /j,  I2  be  two  lines  in  3-space,  let  A,  B  be  two  distinct  points  on  /j,  and  let  C,  D 
be  two  points  on  li-  Let  the  angle  between  the  vertors  AB  and  AC  (resp.  between  AB  and  BD) 
be  a  (resp.  p).  Similarly,  let  the  angle  between  the  vector  CD  and  AC  (resp.  between  CD  and 
BD)  be  -y  (resp.  8).  Then  a>3  implies  that  7>8  (see  Fig.  3.1). 
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Fig.  3.1. 
Proof:  First  note  that  under  these  conditions  we  must  have  C  i^  D  (and  thus  in  particular  7  and 
8  are  well  defined),  for  otherwise  a  would  be  an  interior  angle  in  the  triangle  ABC  which  is 
larger  than  the  exterior  angle  3.  which  is  impossible.   Put  AB  =  u,  i4C  =  x,  BD  =  y.  Then 
CD  =  u  +  y— X.  Since  a>3  we  have 


|x|        ly| 


To  prove  'y>8  we  need  to  show  that 


x(u+y-x)         yju+y-x) 

kl  ly| 

By  the  first  inequality,  it  suffices  to  show  that 
or  that 


|x|  <  |y| 


_  JLl. 


or 


x-y<|xlly| 
which  is  immediate,  since  a  i^  p .  Q.E.D. 

Remark:  The  assertion  in  the  preceding  lemma  can  be  extended  to  the  case  A  =  B  (so  that  a 
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and  P  are  undefined),  li  C  i^  D  then  we  always  have  7  >  8  since  7  is  an  exterior  angle  and  8 
is  another  interior  angle  in  the  triangle  ACD . 

L«nuna  3.3:  The  shortest  path  tt  from  X  to  K  which  passes  through  the  sequence  of  lines 
^1,  .  .  .  ,^„  in  this  order  is  unique.  (Note  that  we  assume  here  that  ^j,  .  .  .  ,4n  are  full  lines, 
or,  alternatively  that  tt  passes  through  interior  points  of  |;i,  .  .  .  ,^„.) 

Proof:  Suppose  that  there  exist  two  shortest  paths  it,  it'  from  X  to  K  through  the  lines 
^1,  .  .  .  ,^„.  Apply  Lemma  3.2  for  each  of  the  (skew)  quadrangles  whose  edges  are  ^,,  ^,+1, 
IT,-,  IT,',  1  =  0,  .  .  .  ,n,  where  tt,-  (resp.  tt,')  is  the  i-th  segment  along  ir  (resp.  it')  (see  Fig. 
3.2). 


Fig.  3.2. 
Lemma  3.1  implies  that  a,+i  =  "y,  and  P,+i  =  8,  for  each  /  =  0,  .  .  .  ,/i  — 1.    The  first 

divergence  of  tt  and  it'  (i.e.  the  first  j  for  which  try  and  IT/'  have  the  same  starting  point  but 

different    endpoints)    forms   a   triangle,    and  we   therefore   have   "y;  >  8/.     Thus,   inductive 

applications  of  Lemma  3.2  imply  that  a,-  >  P<  for  each  i=j+l,  .  .  .  ,n.    In  partiailar, 

a„  >  3„,  which  is  impossible  since  it„  and  ir„'  meet  at  Y,  forming  a  second  triangle  (which 

might  degenerate  to  a  single  segment),  a  contradiction  which  proves  the  lemma.  Q.E.D. 

Lemma  3.3  can  be  strengthened  as  follows:  Call  a  path  it  from  X  to  K  which  passes  through 
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the  lines  4i.  •  •  •  .4n  geodesic  (or  locally  shortest)  if,  for  each  /  =  1,  .  .  .  ,n,  the  path  it 
enters  and  leaves  ^^  at  equal  angles.  It  is  easily  checked  that  every  path  whose  length  is  a  local 
extrerauni  (as  a  function  of  its  points  of  contact  with  ^j,  .  .  .  ,|„)  is  geodesic.  It  is  plain  that 
Lemma  3.3  remains  true  if  one  assumes  lliat  the  paths  in  question  are  only  geodesic.  Hence  we 
have 

Corollary  3.4:  There  exists  a  unique  geodesic  path  ir  from  X  to  Y  which  ipasses  through  a  given 
sequence  of  lines  ^1,  .  .  .  ,^„. 

In  other  words,  the  length  of  the  path  (as  a  funcdon  ol"  the  contact  points)  has  one  global 
minimum,  and  no  other  local  extrema.  Thus,  if  vve  csjntiiiue  to  assume  that  ^j,  .  .  .  .^^  are  full 
lines,  then  we  can  use  two  different  techniques  fov  die  c&lculatiou  of  tlie  contact  points  of  the 
required  shortest  path  with  these  lines.  The  first  tediniquc  ises  approximate  numerical  methods 
for  finding  the  required  minimum.  For  example,  wc  can  irJtially  pass  a  path  tt  through  an 
arbitrary  sequence  of  points,  one  on  each  of  the  given  lines.  Then,  iteratively,  improve  the  path 
by  replacing  each  contact  point  at  which  the  incoming  and  outgoing  angles  are  not  equal  by 
another  point  on  the  same  line  at  which  these  angles  become  equal  (without  changing  the  other 
contact  points).  An  explicit  formula  for  finding  the  new  point  of  contact  can  be  readily  obtained, 
using  elementary  vector  techniques.  (The  problem  involved  here  is,  given  two  points  a  and  b 
outside  a  given  line  /,  to  find  a  point  x  on  /  such  that  both  vectors  x  —  a  and  b  —  x  form  the 
same  angle  with  /.  We  leave  it  to  the  reader  to  verify  that  this  condition  can  be  expressed  by  an 
equation  which  is  quadratic  in  the  coordinates  of  each  of  the  points  a,  b  and  x,  and  whose  overall 
degree  is  4.)  Each  such  iterative  step  shortens  the  length  of  the  path,  and  the  sequence  of  paths 
thus  obtained  will  converge  to  a  path  of  locally  extremal  length,  and  hence  to  the  desired  shortest 
path,  in  virtue  of  Corollary  3.4. 

If  we  wish  to  avoid  numerical  analysis,  and  insist  on  obtaining  a  precise  solution  using  only 
symbolic  calculations,  then  we  can  write  down  a  system  of  n  quartic  equations  m  the  n  positions 
of  the  contact  points  of  the  path  with  the  given  lines,  each  such  equation  corresponding  to  one  of 
the  constraints  given  by  Lemma  3.1,  namely  that  at  each  line  ^  the  incoming  and  outgoing  angles 
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subtended  by  the  shortest  path  be  equal.  This  system  can  then  be  solved  by  elimination  techniques 
(cf.  [Wa]),  leading  to  a  single  polynomial  equation  p(aci)  =  0  in,  say,  the  position  X\  of  the  first 
point  of  contact.  By  CoroOary  3.4,  this  equation  will  have  a  unique  real  solution  which  can  then  be 
rationally  approximated  to  any  desired  degree  of  accuracy. 

Note  however  that  the  resulting  polynomial  p{x)  will  in  general  be  of  degree  which  is  doubly 
exponential  in  n,  because  each  elimination  step  computes  resultants  of  polynomials,  and  thus 
results  in  polynomials  of  one  less  variable  but  of  degree  which  is  roughly  the  square  of  the  degree 
of  the  previous  polynomials,  and  because  n  —  1  such  elimination  steps  are  reqxiired  to  eliminate  all 
but  one  of  the  n  variables  in  question.   This  indicates  that  the  problem  of  computing  the  points  of 
contacts  (and  thus  also  the  length)  of  a  geodesic  path  with  a  sequence  of  line  segments  is  probably 
intractable.  To  be  more  precise,,  suppose  that  we  have  t.wo  candidate  sequences  ^  and  t|  of  lines, 
and  we  wish  to  determine  whether  the  geodesic  path  passing  through  ^  is  longer  than  the  geodesic 
path  passing  through  r\.    This  can  be  solved  precisely  by  standard  methods  involving  symbolic 
calculations  (as  reviewed  e.g.   in  [SS];  see  also  below),  but  these  methods,  which  use  space 
decomposition  techniques  closely  related  to  the  elimination  method  just  noted,  would  also  require 
doubly  exponential  time.    We  have  not  been  able  to  prove  that  this  problem  is  intractable. 
However,  it  is  well  known  that  testing  similar  properties  of  real  roots  of  a  system  of  low  degree 
polynomial  equations  in  n  variables  (or  even  determining  whether  such  roots  exist)  is  NP- 
complete  (this  is  shown  e.g.  by  direct  reduction  from  3-SAT). 

Another  technical  problem  that  should  be  noted  is  that  if  ^j,  .  .  .  ,^„  are  only  segments 
and  not  full  lines,  the  global  minimum  may  be  attained  at  points  lying  outside  those  segments.  In 
this  case  it  is  dear  that  the  shortest  path  it  from  X  to  Y  constrained  to  pass  through  these 
segments  in  order  will  have  to  pass  through  some  endpoints  of  these  segments,  at  which  it  will 
generally  form  unequal  incoming  and  outgoing  angles.  If  this  is  the  case,  and  if  the  enc^xnnts 
through  which  it  must  pass  are  known,  then  the  path-finding  problem  reduces  to  a  collection  of 
subproblems,  each  of  which  calls  for  the  computation  of  the  shortest  path  between  some  pair  of 
points,  which  is  constrained  to  pass  through  a  specified  sequence  of  lines.  The  solution  of  each 
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such  subproblem  can  be  obtained  by  the  methods  outlined  above. 

Since  in  general  the  shortest  path  from  X  to  Y  will  be  a  concatenation  of  subpaths,  each 
connecting  a  pair  of  points  (each  of  which  is  either  X,  Y,  or  a  nonconvex  wall  comer)  and 
constrained  to  pass  through  a  sequence  of  wall  edges,  we  will  '•onsider  the  path-finding  problem  as 
essentially  solved  (up  to  numerical  or  symbolic  calculatians  of  the  sort  discussed  above)  if  we  can 
specify  the  sequence  of  wall  edges  and  comers  through  which  the  desired  shortest  path  must  pass. 
In  this  setting  the  problem  is  reduced  to  a  purely  corobinatoria!  one.  v/hich  we  will  refer  to  as  the 
combinatorial  shortest  path  problem.  It  is  noteworthy  that  this  combinatorial  problem  is  at  least 
solvable  in  finite  time: 

Proposition  3.5:  The  combinatorial  shortest  path  problem  is  solvable  in  doubly  exponential  time 
by  precise,  symbolic  calculations.  It  one  is  allowed  to  use  numerical  analysis  techniques,  then  the 
problem  can  be  solved  in  0{n")  steps,  each  step  consisting  of  finding  a  shortest  path  constrained 
to  pass  through  some  sequence  of  wall  edges  and  comers. 

Proof:  Note  first  that  it  suffices  to  consider  only  a  finite  number  of  possible  sequences  of  wall 
edges  and  comers  through  which  the  shortest  path  from  X  to  y  can  pass,  because  the  shortest 
path  will  not  pass  twice  through  the  same  wall  comer  or  edge.  This  makes  it  plain  that  the  number 
of  such  sequences  that  need  be  considered  is  at  most  0{n").  For  each  pair  of  such  sequences  {, 
T|,  apply  the  procedure  outlined  above  which  uses  symbolic  computations,  to  find  whether  the 
shortest  path  constrained  to  pass  through  the  elements  of  ^  is  shorter  than  the  path  constrained  to 
pass  through  the  elements  of  r\ .  The  sequence  for  which  the  length  of  the  corresponding  path  is 
smallest  is  then  the  solution  to  our  combinatorial  problem.  Note  that  the  basic  step  that  this 
method  employs  is  comparison  between  two  algebraic  numbers,  each  of  which  is  specified  in  terms 
of  the  unique  real  roots  r^,  .  .  .  ,r^  of  some  system  of  Jk  polynomials  in  k  variables,  for  some 
k^n  (these  roots  are  the  points  of  contact  of  one  of  the  shortest  paths  with  the  wall  edges 
through  which  it  is  constrained  to  pass).  Such  a  comparison  can  be  performed  in  precise  terms 
using  Oallins'  cylindrical  algebraic  decomposition  technique  for  analyzing  semi- algebraic  sets 
defined  by  a  system  of  p  polynomial  equalities  and  inequalities  in  k  variables,  having  maximum 
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degree  m  ([Co];  see  also  [SS]).  Collins'  technique  runs  in  0((mp)  )  time.  In  our  case  each 
polynomial  equation  is  quartic,  andp,  Jk  ^  n.  It  follows  that  the  overall  cost  of  the  combinatorial 
shortest  path  problem  is  doubly  exponential  in  n.  When  numerical  methods  are  allowed  in  the 
evaluation  of  each  of  the  shortest  paths  constrained  to  pass  through  some  sequence  of  wall  edges 
and  comers,  the  problem  can  plainly  bie  solved  in  0(^n")  such  numerical  evaluations.  Q.E.D. 

It  is  currently  an  open  problem  whether  faster  procedures  than  the  straightforward  one  just 
sketched  exist  for  solving  the  combinatorial  shortest  path  problem. 

4.  Shortest  Paths  Along  a  Convex  Polyhedron 

In  this  section  we  analyze  the  problem  of  calculating  the  shortest  path  between  two  points 
along  the  surface  of  a  convex  polyhedron  In  3-space.  This  special  case  is  favorable  because  of 
various  properties  of  geodesic  paths  along  a  convex  polyhedron.  These  properties  will  be  analyzed 
in  this  section;  an  algorithm  for  the  calculation  of  such  shortest  paths  will  be  presented  in  the 
following  section. 

Let  AT  be  a  given  convex  polyhedron,  and  let  S  denote  its  boundary.  Let  X  and  Y  be  two 
points  on  S.  The  problem  that  we  consider  is  to  calculate  the  shortest  path  from  X  to  y 
constrained  to  lie  along  S.  It  will  be  more  convenient  to  consider  a  somewhat  more  general 
problem,  namely  -  given  a  point  X  on  S,  we  wish  to  pre-process  K  so  that  later,  for  any  desired 
destination  point  y  on  S,  the  shortest  path  from  X  to  K  can  be  easily  and  quickly  calculated.  To 
simplify  the  foregoing  analysis  we  will  assume,  without  real  loss  of  generality,  that  the 
representation  of  ^  is  nondegenerate,  in  the  sense  that  no  two  faces  oi  K  art  coplanar. 
(Otherwise,  we  can  repeatedly  combine  pairs  of  adjacent  co-planar  faces  into  single  faces,  until  the 
above  property  holds  for  K.)  Let  n  be  the  number  of  vertices  of  K,  so  that  K  has  also  0(n) 
edges  and  faces. 

Definition:  (a)  A  point  Z  €  5  is  called  a  ridge  pcnnt  if  there  exist  at  least  two  shortest  paths 
from  X  XoZ  along  S  (cf.  Fig.  4.1.)  We  denote  by  i?  the  set  of  all  ridge  points  in  S. 

(b)  For  each  point  Z  ^  5  —  i? ,  let  Tr(Z)  denote  the  unique  shortest  path  from  X  XoZ  along  S. 
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Fig.  4.1.  Ridge  points  along  a  convex  polyhedron 
We  will  first  cx)nsider  the  following  subproblem:  Let  Z  €  5,  and  suppose  that  the  sequence 

5  =  (§1,  .  .  .  ,4„)  of  edges  of  K  through  which  Ti{Zi)  passes  has  already  been  found  (we  will 

shortly  see  in  Lemma  4.1  that  17(2)  cannot  cross  a  vertex  of  K,  so  that  it  meets  each  of  the  edges 

^l»  •  •  -  »4rt  in  an  interior  point  of  that  edge).  Then  we  wish  to  Qnd  the  points  of  contact 

between  tt(Z)  and  each  of  the  edges  ^,.  In  the  general  3-dimensional  case  treated  in  the  previous 

section,  the  solution  of  this  subproblem  had  been  rather  complicated,  and  involved  solution  of 

high-degree  algebraic  equations,  mainly  because  any  two  adjacent  edges  4,_i  and  ^,  in  ^  could  be 

skew  to  one  another.    However,  in  the  special  case  which  concerns  us  here  this  cannot  happen, 

thereby  making  the  problem  immediately  solvable  as  follows. 

Let  IT  consist  of  the  segments  tto,  .  .  .  ,1T„.  Recall  that  Lemma  3.1  implies  that  for  each 
J  =  l,  .  .  .  ,n,  the  angles  that  tt,_j  and  ttj  subtend  at  ^,  are  equal.  This  suggests  the  following 
simple  algorithm  for  the  calculation  of  the  points  of  contact.  Since  we  know  the  sequence  of  edges 
through  which  Tr(Z)  passes,  we  also  know  the  corresponding  sequence /q/i,  .  .  .  ,f„  of  faces 
of  K  through  which  'ir(Z)  passes,  where  the  face  f(  contains  the  two  edges  ^/  and  5,+i,  for 
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i  =  l,  .  .  .  ,n  —  l,  and  where /q  contains  X  and  ^1  and /„  contains  ^„  and  Z. 

We  then  unfold  the  collection  of  faces  /q,  ...,/„  so  as  to  make  them  all  lie  in  the  same  plane 
L.  This  is  done  iteratively.  That  is,  initially  we  place /q  in  L,  letting  X  coincide  with  the  origin. 
Suppose  that  we  have  already  unfolded  and  placed  in  L  all  faces  up  to/,_i.   We  then  unfold/, 
about  ^j  until  it  becomes  co-planar  with  /,_i  (but  lies  on  the  other  side  of  ^,).    In  practice,  we 
compute  for  each  face  /,  the  displacement  a,-  and  orientation  6,  defining  its  position  in  L  relative 
to  some  standard  and  fixed  plane  representation  of  this  face.  We  can  then  compute  from  a„  and 
6„  the  position  of  Z  in  L.  The  required  path  tt(Z),  unfolded  to  L,  is  then  simply  the  straight 
segment  XZ.  The  points  of  intersection  of  this  segment  with  the  unfolded  edges  ^i,  .  .  .  ,  ^„  are 
then  readily  determined,  and  can  be  easily  traniformed  back  to  the  original  polyhedron.    For 
further  reference,  let  us  call  this  process  as  the  planar  urrfolding  of  K  relative  to  ^i,  .  .  .  ,^„;v/c 
also  refer  to  the  pair  {a„,Q„)  as  the  position  of /^  in  that  planar  unfolding.   (See  also  Alexandrov 
[Al]  for  an  analysis  of  the  unfolded  planar  structure  of  AT;  the  above  observations  have  also  been 
made  by  Franklin  et  al.  [FAV],  [FA]  although  they  have  not  developed  them  into  a  polynomial- 
time  algorithm). 

Hence,  as  in  the  general  3-D  case,  the  main  problem  which  needs  to  be  solved  is  that  of 
calculating  the  sequence  ^j,  .  .  .  ,^„  of  edges  through  which  the  path  ir(y)  passes. 

To  this  end,  we  will  partition  S  into  at  most  n  vertex-free  connected  regions,  caHedpeeb,  such 
that  the  interiors  of  these  regions  do  not  contain  any  ridge  point,  and  such  that  for  each  such 
region  p,  the  path  'It(Z)  to  any  Z  ^  p  h  wholly  contained  in  p.  Since  p  is  vertex-free,  the 
sequence  of  edges  through  which  Tt(Z)  passes  will  be  easy  to  calculate,  as  will  be  shown  below. 

To  obtain  this  partitioning,  we  begin  with  aiudysis  of  several  properties  of  ridge  points. 

Lemma  4.1:  A  shortest  path  'n'(Z)  caimot  pass  through  a  vertex  of  K. 

Proof:  (We  are  indebted  to  R.  Pollack  for  suggesting  this  simplified  proof.)  Suppose  the  contrary, 
and  let  6''  be  a  vertex  of  K  lying  on  a  shortest  path  'Tr(Z)  from  X  to  some  Z  ^  S.  Suppose  first 
that  K  is  incident  to  exactly  three  faces  of  K.  Let  A  and  B  be  two  points  on  'n{Z)  lying  on  the 
two  straight  subsegraents  of  tt(Z)  adjacent  to  U,  with  A  lying  before  U  and  B  after  U  along  that 
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path.  Let  f/^.f^  be  the  faces  containing  A  and  B  respectively,  and  let  /  be  the  third  face  of  K 
containing  V .  Instead  of  the  planar  unfolding  of  iC  in  which  'Tr(Z)  is  a  straight  segment  (call  this 
the  "straight"  unfolding  of  17(2)),  we  construct  another  unfolding  as  foUows.  First  unfold  faces 
of  K  into  the  plane  as  in  the  planar  unfolding  of  'ir(Z),  until  U  is  reached.  Then,  instead  of 
unfolding /5  past/^  (as  is  done  in  the  straight  unfolding  of  17(2)),  unfold  first/  past/^  along 
their  common  edge,  and  then  unfold  /^  past  /  along  th^  common  edge;  then  continue  to  unfold 
as  in  the  remainder  of  the  straight  unfolding  of  it(Z)  (see  Fig.  4.2). 
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Fig.  4.2.  Proof  of  Lemma  4.1. 
In  this  new  unfolding  the  path  ir(Z)  appears  as  a  broken  line  at  U ,  and  it  is  easily  seen  that 

'it(Z)  can  be  shortcut  near  C/,  yielding  a  shorter  path  which  is  also  contained  within  this  new 

unfolding.  This  contradiction  establishes  the  lemma.  Essentially  the  same  argument  also  applies  in 

case  U  is  incident  to  more  than  three  faces.   Q.E.D. 

Lemma  4.2:  A  shortest  path  Tr(Z)  cannot  pass  through  a  ridge  point. 

Proof:  Suppose  that  tt(Z)  does  pass  through  a  ridge  point  ^ .  Then  the  initial  portion  tti  of 
'It(Z)  up  to  W^  is  one  of  several  shortest  paths  from  X  to  W.  Let  ir2  be  another  »uch  shortest 
path,  and  without  loss  of  generality  assume  that  ttj  and  112  are  transver»al  to  one  another  at  TV 
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(this  will  be  the  case  if  we  choose  W  to  be  the  first  ridge  point  along  it(Z),  which  is  always  well 
defined,  because  the  set  or  ridge  points  along  tt(Z)  is  dosed,  as  can  be  easily  verified).  Let  tt' 
denote  the  path  obtained  by  replacing  tTj  by  it 2  in  tt(Z);  note  that  ir'  is  also  a  shortest  path  to 
Z.  However,  if  W  is  interior  to  some  face  of  K,  then  it'  cannot  be  a  shortest  path  to  W,  since  it 
bends  at  an  interior  point  of  a  face  of  K.  On  the  other  hand,  if  W  lies  on  an  edge  of  K,  then  it' 
cannot  be  a  shortest  path  to  W  since  it  forms  unequal  angles  with  the  edge  containing  W  (note 
that  W  cannot  be  a  vertex  by  Lemma  4.1).  Q.E.D. 

Lemma  4.3:  The  set  R  of  ridge  points  is  the  union  of  finitely  many  straight  segments. 

Proof:  With  any  point  Z  of  /?  we  can  associate  the  face  of  K  containing  Z,  and  the  two  sequences 
of  edges  of  K  through  which  the  two  shortest  paths  from  X  to  Z  pass.  Since  there  are  only  finitely 
many  values  that  each  of  these  three  parameters  can  assume  (because  a  shortest  path  to  a  point 
cannot  pass  through  an  edge  of  K  more  than  once),  it  suffices  to  show  that  the  locus  of  all  ridge 
points  Z  lying  on  a  fixed  face  of  K,  for  which  the  two  shortest  paths  from  X  to  Z  pass  through 
two  fixed  sequences  of  edges  of  K,  is  a  straight  segment. 

Let  therefore/  be  a  fixed  face  of  K,  and  let  4  =  (^1.  •  •  •  .^n).  "H  =  ('ni»  •  •  •  •'Hm)  ^  ^° 
fixed  sequences  of  edges  of  K,  such  that  any  two  adjacent  edges  in  either  sequence  lie  on  a 
common  face,  and  such  that  ^i  and  1)1  bound  the  face  containing  X,  while  ^„  and  T|^  bound/. 
Assume  that  there  exists  at  least  one  point  Z  ^  R  having  these  values  as  its  associated 
parameters.  Let  (a^,9^),  {a^,Qy^)  be  the  positions  of  /  in  the  planar  unfoldings  of  K  relative  to 
the  two  sequences  ^  and  T|  respectively,  where  a^  and  a_  both  give  the  position  of  Z  in  the 
corresponding  planar  unfolding.  Let  W  ^  R  be  another  point  having  the  same  parameters  as  Z, 
and  write  ZW  =  w  in  the  standard  Cartesian  representation  of  /.  Then  we  have  (where  Re 
(resp.  R^)  denotes  the  rotation  of  the  plane  by  the  angle  6c  (resp.  6._)): 

which  states  the  equality  of  the  lengths  of  the  two  shortest  paths  from  X  to  W,  constrained  to 
pass  respectively  through  the  sequences  ^  and  i\  of  edges  of  K.  Squaring  out  the  above  equation, 
and  using  the  fact  that  |^^w|  =  [/Jt^wI,  and  that  |a^|  =  |a.^|,  we  obtain 
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which  is  the  equation  of  a  straight  line  passing  through  Z.  Q.E.D. 

Remarks:  (1)  The  preceding  equation  is  nondegencrate  provided  that  fl^c  ^  ^tj^ti-  But  if  these 
two  quantities  were  equal,  then  the  two  corresponding  shortest  paths  to  W  would  be  such  that 
their  terminal  segments  along  /  coincide.  This  however  implies  that  these  two  paths  pass  through 
another  ridge  point  (e.g.  any  interior  point  of  these  coincident  segments),  which  is  impossible,  by 
Lemma  4.2. 

(2)  The  equation  defining  the  ridge  segment  given  above  can  be  rewritten  as 

(a^g  -  a^^)'W  =  0 

and  thus  be  given  the  foUowing  geometric  interpretation: 

Perform  the  two  planar  unfoldings  relative  to  ^  and  T|.  First  rotate  each  unfolding  in  clockwise 
direction  by  the  respective  angle  9c,  6_  (note  that  these  rotations  cause  the  final  face/  in  each  of 
the  unfoldings  to  have  the  same  orientation  as  the  standard  representation  of/).  Next  translate 
(without  rotating)  the  two  resulting  copies  of  /  so  that  they  both  coincide  with  the  standard 
representation  of  /.  Note  that  this  will  have  moved  the  starting  point  X  to  the  planar  positions 
X^  =  fl^i  and  X^  —  a-ri^Ti  respectively.  Then  the  required  ridge  segment  (in  the  standard 
representation  of/)  is  contained  in  the  perpendicular  bisector  to  the  segment  X^^  (sec  Fig.  4.3). 
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Fig.  4.3.  Qanstniction  of  ridge  segments 
(3)  This  construction  of  ridge  segments  makes  it  dear  that  if  Z  lies  on  a  ridge  segment  e  defined 
by  two  edge  sequences  ^  and  t\  and  Z  is  not  a  vertex  of  K,  then  for  points  W  lying  on  the  line 
containing  e  in  a  sufficiently  small  neighborhood  of  Z  the  following  properties  hold:  There  exist 
two  geodesic  (but  not  necessarily  shortest)  paths  to  W  having  equal  lengths  and  passing 
respectively  through  the  sequences  ^  and  x\\  moreover  the  starting  orientation  of  each  of  these 
paths  is  a  monotone  and  continuous  function  of  W.  (Note  that  these  geodesies  will  be  shortest 
paths  for  W  lying  on  at  least  one  side  of  Z.) 

To  continue  our  analysis  of  the  structure  of  /? ,  we  first  introduce  the  following  notation. 

Drflnltion:  (a)  For  each  planar  orientation  9  define  a  polygonal  path  p  ■■  p(0)  from  X  along  S 
as  follows:  p  starts  at  X  in  the  direction  of  9  (relative  to  the  standard  representation  of  the  ixcf^ 
containing  X).  Whenever  p  reaches  an  edge  tf  of  i^  it  bends  over  it  to  the  adjacent  face  so  that 
the  two  segments  of  p  adjacent  to  e  form  with  it  equal  angles  (as  in  Lemma  3.1).  p  will 
terminate  as  soons  as  it  reaches  a  vertex  of  K. 
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(b)  For  each  orientation  6  and  each  Z  €  p(6)  we  define  p(9^e)  to  be  the  initia]  portion  of 
p(6)  between  X  and  Z. 

With  some  exceptions  noted  below,  we  will  consider  the  path  p(9)  to  terminate  after  the 
first  time  it  reaches  either  a  vertex  of  /iT  or  a  ridge  point.  (Note  that  this  must  occur  for  each  6, 
because  otherwise  either  the  length  of  p  would  inaease  without  bound,  while  p  still  being  the 
shortest  path  to  any  of  its  points,  which  is  plainly  impossible,  or  else  p  would  reach  X  again,  in 
which  case  it  is  clear  thatp  contains  a  ridge  point.)  Ltt  r(0)  denote  the  endpoint  of  p(9). 

Lemma  4.4:  (a)  Let  Z„  be  a  sequence  of  jx»ints  on  S  converging  to  some  Z  ^  5  as  n  -•  <».  For 
each  n  ^  1  let  IT  „  be  a  shortest  path  from.  X  to  Z„,  aad  suppose  that  the  paths  ir^  converge  in 
the  Hausdorff  topology  of  sets  to  a  path  it.  Then  ir  is  a  shortest  path  to  Z. 

(b)  The  function  r(9)  is  continuous,  ;    -    ■ 

Proof:  (a)  Suppose  the  contrary,  and  let  ir'  be  a  path  from  X  to  Z  which  is  shorter  than  it.  Since 
the  length  of  ir  is  the  limit  of  the  lengths  of  the  paths  ir^,  it  follows  that  if  n  is  sufficiently  large, 
by  appending  to  ir'  a  short  path  connecting  Z  to  Z„,  we  can  obtain  a  path  to  Z„  which  is  shorter 
than  7r„,  a  contradiction. 

(b)  Let  9„  -  9,  and  suppose  that  Z„  =  r(9„)  converges  to  some  point  Z,  and  that  for  each  n 
the  point  r(9„)  is  a  ridge  point  (if  the  latter  property  cannot  be  achieved,  then  from  a  certain  n 
on,  r(Q„)  is  a  vertex,  which  implies  that  the  sequence  9„  is  constant  too,  in  which  case  there  is 
nothing  to  prove).  It  follows  that  for  each  n  there  exists  another  orientation  <{>„  such  that 
'■(^n)  ~  ''{^n)  ~  2„.  Passing  to  subsequences  if  necessary,  we  can  assume  that  4>„  converge  to 
some  (j>,  and  that  the  paths  p(9„,Z„)  converge  to  p{B,Z)  and  similarly  the  paths  p(4>„^„) 
converge  to  p(4>,Z).  Let  ^'*,  t)"  be  the  two  sequences  of  edges  of  K  through  which  the  two 
paths  p{Q„^„)  and  p(^^,Z„)  pass.  Passing  again  to  a  subsequence  if  necessary,  we  can  assume 
that  the  sequences  ^"  and  t)"  are  both  constant  (necessarily  distinct  from  one  another).  Two 
cases  can  arise: 

(i)  If  9  ?^  4>  then  by  (a)  the  two  distinct  paths  p(Q,Z)  and  p(<}>yZ)  are  shortest  paths  to  Z,  so 
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that  Z  is  a  ridge  point,  and  r(6)  =  Z,  because  p(6yZ)  does  not  pass  through  a  vertex  of  ^  or  a 
ridge  point  (other  than  Z). 

(ii)  If  6  =  <{)  then  Z  must  be  a  vertex  of  K.  Indeed,  if  this  were  not  the  case,  then  (by  Lemma 
4.1)  there  would  exist  a  sufficiently  small  neighborhood  U  of  p(0,Z)  which  contains  no  vertex  of 
K.  But  then  for  n  sufficiently  large  the  two  paths  p(9„^„),  p(^„^„)  would  be  wholly 
contained  in  U,  and  thus  will  cross  the  same  sequence  of  edges  of  K,  contrary  to  assumption.  The 
same  argument  as  in  (i)  above  now  implies  that  Z  =  r(9). 

This  proves  part  (b)  of  the  lemma.  Q.E.D. 

Corollary:  The  set  /?  *  consistiiig  of  all  vertices  of  K  and  ridge  points  is  a  dosed  connected  set. 

Proof:  It  suffices  to  show  that  the  map  r  defined  above  is  onto  R*,  tor  then  R *  will  be  the 
continuouis  image  of  the  unit  circle  of  starting  orientations  6.  To  show  that  r  is  onto,  let  Z  €  R*, 
and  let  'It(Z)  be  one  of  the  shortest  paths  to  Z.  Then,  arguing  as  in  the  preceding  proof, 
Z  =  '■(9),  where  6  is  the  starting  orientation  of  ir(Z)  (as  follows  from  Lenunas  4.1  and  4.2  and 
from  the  definition  of  r).  Q.E.D. 

In  view  of  Lemma  4.3,  the  set  /?  *  can  be  regarded  as  a  graph  whose  edges  are  (portions  of) 
the  straight  segments  yielded  by  the  proof  of  the  preceding  Lemma.  The  vertices  of  tliis  graph  are 
either  vertices  of  K,  or  points  at  which  such  a  segment  intersects  an  edge  of  K,  or  points  at  which 
two  sudi  segments  intersect.  The  degree  of  each  vertex  u  of  R*  is  defined  to  be  the  number  of 
edges  of  /?*  incident  to  u.  More  information  on  the  structure  of  the  vertices  of  /?  is  obtained 
from  the  following  lemmas. 

Lenmia  4.5:  (a)  Each  vertex  oi R*  having  degree  1  is  a  vertex  of  K. 

(b)  Eadi  vertex  of  /?  *  having  degree  2  is  an  intersection  of  R  with  the  interior  of  some  edge  of 
K. 

(c)  R*  does  not  contain  any  closed  path. 

Proof:  (a)  Suppose  that  Z  is  a  vertex  oi  R*  of  degree  1  which  is  not  a  vertex  of  K.  Then  Z  must 
be  a  ridge  point.   Suppose  first  that  there  are  exactly  two  shortest  paths  it  and  ir'  reaching  Z.  By 
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Lemma  4.3  Z  lies  on  a  segment  e  which  is  determined  by  the  two  seqiiences  t>  ^  o^  edges  of  K 
through  which  ir  and  it'  respectively  pass.  If  Z  is  an  endpoint  of  e  (and  of  no  other  ridge 
segment)  then  for  each  point  W  on  the  line  containing  e  lying  near  Z  on  the  other  side  of  e  there 
must  exist  a  unique  shortest  path  it(IiV)  to  W.  However,  this  path  passes  through  a  sequence  { 
of  edges  of  K  which  is  distinct  from  both  ^  and  t).  For  suppose  that  'ir(W)  passes  through  the 
sequence  5;  then  by  definition  of  e  there  would  east  another  path  from  X  to  W  having  the  same 
length  as  'It(W)  and  passing  through  the  sequence  t|  of  edges  (cf.  Remark  (3)  foUowing  Lemma 
4.3).  Hence,  letting  W  -  Z  along  this  line,  we  wo'Jd  obtain  a  third  shortest  path  to  Z  passing 
through  the  sequence  5,  contrary  to  assumption.  This  argument  can  be  generalized  to  the  case  in 
which  more  than  two  shortest  paths  reach  Z,  thus  proving  (a). 

(b)  Suppose  the  contrary,  and  let  Z  €  i?  be  an  interior  point  of  some  face  /  of  K,  which  is  a 
common  endpoint  of  exactly  two  ridge  segments  e\  and  e^-  Let  ^  ,  T]  be  two  distinct  pairs  of 
sequences  of  edges  of  K  which  define  c,,  J  =  l,2.  Suppose  for  definiteness  that  ^  ^  4  .  Then 
there  are  at  least  three  shortest  paths  reaching  Z  and  passing  respectively  through  the  sequences 
^\  ^  and  T^  (see  Fig.  4.4).  Consider  the  perpendicular  bisector  /  to  the  segment  XciXcZ  as  in 
the  ridge  construction  procedure  described  above.  /  is  distinct  from  the  lines  containing  e\  and 
e2i  it  passes  through  Z,  and  contains  no  ridge  point  in  a  small  neighborhood  of  Z,  except  Z  itself. 
Thus  for  each  W  ^  /  near  Z  there  exists  a  unique  shortest  path  'it(W)  passing  through  a 
sequence  C(W^)  of  edges  of  K  which  is  distinct  from  ^',  r} ,  I  =  1,2,  as  can  be  seen  from  Remark 
(3)  following  Lemma  4.3.  But  then,  using  arguments  similar  to  those  used  in  the  preceding 
jjaragraph,  we  obtain  a  contradiction  which  concludes  the  proof  of  (b). 
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Fig.  4.4.  Illustrating  the  proof  that  ridge  segments 
cannot  "bend"  at  a  point  interior  to  a  face. 

(c)  SuRX)se  the  contrary,  and  let  C  be  a  simple  dosed  path  in  /? *.  By  the  Jordan  cxirve  theorem, 
C  divides  5  into  two  disjoint  regions  S^,  S2,  so  that  one  of  them,  say  Si  contains  X.  Let  Z  be  an 
arbitrary  point  in  82-  Then  the  shortest  path  1^(2)  from  X  to  Z  will  have  to  mtersect  C,  which 
contradicts  either  Lemma  4.1  or  Lemma  4.2.  Q.E.D. 

CoroOary:  R*  is  a  tree  having  (some  of)  the  vertices  of  Jf  as  leaves. 
Proof:  Immediate  by  the  preceding  lemma. 

Suppose  that  we  have  managed  to  construct  R  * .  We  can  then  use  it  for  calctilating  shortest 
paths  from  X  to  arbitrary  points  Y  ^  S  as  follows: 

Etefine  Q  to  be  the  union  of  R  with  the  shortest  paths  from  X  to  every  vertex  of  K.  Q  partitions 
S  into  disjoint  connected  regions  (which  we  call  peels),  the  interiors  of  which  do  not  contain  any 
vertex  or  ridge  point.  Let  p  be  one  of  these  peels.  Using  arguments  similar  to  the  preceding  ones, 
we  can  show  that  the  (unique)  shortest  path  'n'(Z)  to  any  interior  point  Z  of  p  is  wholly  contaiived 
in    p.      Also,     there     exists     an     open     interval    /  =  /(p)     of     orientations     Mich     that 
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p  =     U    p(6,r(6)),  and  such  that  for  each  9  €  /,  the  endpoint  r(6)  docs  not  reach  a  vertex. 

e  €  / 

Moreover,  the  two  end  orientations  8j  and  82  of  /  are  such  that  bothp(6i)  and/7(82)  terminate 
at  a  vertex  of  K  (otherwise  /  and  p  could  be  increased  with  all  the  other  properties  of  p  still 
valid).  Finally,  the  intervals  /(/>)  are  pairwise  disjoint,  and  thie  union  of  their  dosures  cover  the 
whole  angular  space  [0,211]. 

Let  /„  denote  the  set  of  edges  of  K  which  intersect  p.  It  follovv'S  that  the  edges  /_  can  be  ordered 
by  their  adjacency  in  p.  That  is,  we  say  that  e  precede  e'  in  j'  if  there  exists  a  shortest  path  in  p 
passing  first  through  e  and  later  through  ^' .  To  see  that  this  order  is  well  defined,  and  to  gain 
more  insight  into  the  structure  of  peels,  we  have  the  following  lemmas: 

Lemma  4.6:  Suppose  that  a  peel  p  is  unfolded  into  the  plane  by  u-ofolding  each  of  the  geodesic 
paths  p(6)  in  p  into  the  ray  from  X  at  orientation  6.  Then  the  resulting  image  of  p  is  convex. 

Proof:  Let  f  be  the  unfolded  image  oi p.  P  is  plainly  contained  in  the  angular  sector  consisting  of 
rays  at  orientations  6  d  (9i  ,82)  and  is  star  shaped  with  respect  to  X.  Hence  if  P  is  not  convex 
there  must  exist  a  concave  comer  Z  on  its  boundary  but  not  on  either  of  the  two  rays  9  =  6^  or 
8  =  82.  Hence  Z  is  a  ridge  point.  Let  the  two  (ridge)  edges  meeting  at  Z  be  « j  and  «2.  ^nd  let 
/  l)e  the  line  containing  e^.  Then  points  W  ^  I  near  Z  but  on  the  other  side  of  Cj  will  be  inside 
P,  so  that  the  shortest  path  to  such  a  W  is  unique,  is  contained  in  p,  and  passes  through  the  same 
sequence  of  edges  as  the  path  which  reaches  Z  from  within  p.  But  then  Remark  (3)  following  the 
construction  of  ridge  segments  given  above  implies  that  there  also  exists  another  geodesic  path  to 
W  whose  length  is  equal  to  the  length  of  the  geodesic  path  to  W  through  p .  This  contradiction 
proves  the  assertion.  Q.E.D. 

Lemma  4.7:  Let  p  be  a  peel  and  let  /  be  a  face  of  K  such  that  pH/  ^  0.  Then  there  exists  a 
unique  sequence  ^  of  edges  of  K  such  that  for  each  Z  €  pH/  the  shortest  path  'It(Z)  to  Z 
passes  through  the  sequence  |. 

Proof:  Suppose  the  contrary,  and  let  Z^,  Z2  be  two  points  in  pfl/  for  which  the  sequences  {,  "H 
of  edges  of  K  through  which  the  paths  tt(Zi),  'ir(Z2)  respectively  pass  are  distinct.   By  Lemma 
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4.6  the  segment  Z1Z2  is  pontained  in  pfl/,  and  by  continuity  there  would  have  to  exist  a  ridge 
point  Z  ^  Z1Z2  (hence  inp),  a  contradiction  which  proves  our  claim.  Q.E.D. 

Hence  each  peel  p  defines  at  most  0(/l)  distinct  sequences  of  edges  of  K,  one  for  each  face 
of  K  whir'a  p  intersects,  such  that  the  shortest  path  to  any  Z  ^  p  passes  through  one  of  these 
sequences.  These  sequences  can  be  arranged  in  an  auxiliary  tree  associated  with  p  so  that  ^  is  an 
ancestor  of  i]  in  this  tree  if  and  only,  if  ^  is  a  prefix  of  t].  Given  a  peel  p,  it  is  straightforward  to 
compute  all  these  sequences,  and  we  omit  details  of  this  easy  construction. 

These  observations  lead  to  the  foUowing 

Proposition  4.8:  (a)  There  are  n  peels. 

(b)  There  are  O(n^)  edges  xnR*.'   -  -^  ■■       -■ 

Proof:  (a)  Immediate  from  the  preceding  considerations. 

(b)  We  will  show  that  each  face  of  K  contains  0(n)  segments  of  /?*.  Let  /  be  a  face  of  K,  and 
let  ^^  ...  ,^'  denote  the  collection  of  all  sequences  of  edges  of  K  traversed  by  shortest  paths 
from  X  to  points  on/.  Each  such  sequence  corresponds  to  a  unique  peel,  so  that,  by  Lemma  4.7, 
t  ■^  n.  For  each  such  sequence  ^'  let  (a,/, 9c/)  be  the  parameters  describing  the  position  of  /  in 
the  planar  unfolding  corresponding  to  the  sequence  ^  ,  and  let  Xct  =  dcJi-^    denote  the  planar 

position  of  the  point  X  when  this  unfolding  is  moved  so  that  /  coincides  with  its  standard  planar 
representation. 

In  a  manner  quite  similar  to  the  analysis  of  Voronoi  diagrams  [Sh]  we  can  define  the  dual  graph 
of  R*nf  to  consist  of  the  points  Xci  as  its  nodes,  such  that  Xd  and Xrj  are  connected  by  an  edge 
if  the  two  sequences  ^'  and  ^  define  an  edge  of  R*nf.  As  in  [Sh],  one  can  show  that  this  dual 
graph  is  planar  by  embedding  it  into  the  plane  as  follows.  Map  each  node  Xr'  to  its  plane 
position,  and  map  each  edge  {Xd,  Xd)  to  the  union  of  two  segments  connecting  respectively  Xd, 
Xd  to  a  point  on  the  common  edge  of  R  Of  which  these  two  sequences  define.  Since  this  graph 
is  planar,  and  since  Lemma  4.5(b)  implies  that  it  has  no  multiple  edges,  it  foUows  that  it  has  at 
most  0{n)  edges,  thus  proving  (b).  (A  similar  connection  between  the  partitioning  of  5  by  the  set 


-28- 

of  ridge  points  and  Voronoi  diagrams  is  also  noted  in  [FAV].)  Q.E.D. 

Remark:  The  slicing  of  the  surface  S  of  K  into  peels  has  the  following  geometric  implication.  If 
we  apply  the  planar  unfolding  procedure  to  all  the  peels,  from  the  same  planar  position  of  X,  we 
obtain  a  planar  layout  which  we  denote  as  U{K)  which  has  the  following  properties: 

(1)  No  two  unfolded  peels  overlap  in  the  plane  (except  at  points  lying  on  an  unfolded  geodesic 
path  connecting  X  to  a  vertex  of  K).  Indeed,  if  Z^  is  a  point  in  one  peel  and  Z2  is  a  point  in 
another,  then  the  segments  XZj  and  XZ2  have  different  orientations,  unless  the  exceptional 
condition  noted  above  holds. 

(2)  U{K)  is  star  shaped  with  respect  to  X. 

(3)  If  V(K)  is  another  planar  unfolding  of  the  whole  surface  of  K  having  properties  (1)  and  (2), 
then  the  smallest  disc  about  X  containing  U{K)  has  radius  smaller  than  or  equal  to  that  of  the 
corresponding  such  disc  containing  V(K). 

The  partitioning  of  S  into  peels  by  Q  enables  us  to  find  the  shortest  path  from  X  to  any 
K  €  5  in  the  following  simple  manner: 

(a)  Find  the  peel  p  containing  Y.  For  simplicity  assume  that  7  is  an  interior  point  of  p. 

(b)  Fmd  the  sequence  ^  of  edges  in  !„  through  which  17(7)  passes. 

(c)  Apply  the  planar  unfolding  procedure  relative  to  ^,  as  described  above,  to  obtain  the  required 
shortest  path  tt(Y). 

ODnceming  the  complexity  of  the  path-finding  procedure  just  outlined,  assume  that  the  face 
f  oi  K  containing  Y  is  specified  too.  Each  such  /  is  partitioned  by  Q  into  at  most  0{n) 
subregions.  It  is  then  straightforward  to  locate  in  time  0{n)  the  subregion  of/  containing  Y,  and 
thus  the  peel  p  containing  Y  (for  example,  pass  a  straight  line  through  Y  and  compute  its 
intersections  with  QDf,  from  which  the  region  containing  K  can  be  readily  calculated),  dveap, 
we  obtain  ^  in  0(n)  time  from  the  pre-calculated  associated  tree,  and  the  planar  unfolding 
procedure  can  then  be  applied  to  ^  also  in  time  0(jl). 
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CoroUary:  After  pre-processing,  the  shortest  path  from  a  fixed  point  X  to  any  point  Y  ^  S  can  be 
computed  in  time  0{n),  using  a  data-structure  whose  size  is  0(rr). 

The  main  problem  that  we  still  face  is  that  of  computing  Q.  This  will  be  done  in  the 
following  section. 

5.  The  Peels  Construction  Algorithm. 

In  this  section  we  present  an  algorithm  for  partitioning  the  surface  of  K  into  peels,  which 
runs  in  time  0{n  log  n).  The  algorithm  constructs  a  tree  T,  called  the  slice  tree  as  follows.  Let 
us  define  a  slice  a  in  terms  of  two  "starting  OTientatiom"  61O2  and  a  "terminal  face"/,  which 
are  assumed  to  have  the  following  properties: 

(i)      Each  of  p(9i)  and p (62)  reaches  a  vertex  of  K  either  before  readiing/  or  on/  itself. 

(ii)  For  each  6  €  (01,62)  ^^  P^'^  P(6)  reaches  /  and  does  not  meet  a  vertex  before  its  first 
exit  from/. 

(iii)  For  each  6  ^  (61,62)  let  Ze  €  p(6)  be  the  point  at  which  this  path  leaves  /  (for  the  first 
time);  then  all  the  points  Zq  lie  on  the  same  edge  of  /  (which  we  call  the  terminal  edge  of 
cr). 

The  corresponding  slice,  denoted  as  O'(6i,62v0>  ^  defined  to  be  the  union  of  all  paths 
P{Q^q),  for  6  €  (61,62).  Note  that  the  set  of  points  at  which  geodesic  paths  in  o-(6i,62/) 
enter  /  is  also  a  subsegment  of  some  edge  of/.  Note  also  that  all  the  geodesic  paths  p(9^e)  in 
this  slice  pass  through  the  same  sequence  of  edges  of  K,  and  that  if  ct  is  unfolded  into  the  plane 
along  this  sequence  of  edges,  its  image  is  a  triangle  bounded  by  the  two  rays  emerging  from  the 
origin  (X)  at  orientations  9i  and  62  and  by  the  terminal  edge  of  cr  (cf.  Fig.  5.1). 
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Fig.  5.1.  Planar  layout  of  a  slice. 
A  slice  is  roughly  meant  to  correspond  to  some  portion  of  a  peel,  which  is  'cacodcd' 

implicitly  by  the  range  of  starting  orientations  and  the  terminal  face/.  The  correspondence  is  such 

that  for  each  'true'  slice  a  =  cr(9i,02x/)  the  sequence  of  edges  that  it  meets  is  one  of  the 

sequences  associated  with  some  peel  p  (in  the  manner  described  in  the  preceding  section),  and  the 

range  (61,62)  is  contained  in  the  corresponding  range  I(p)-  Note  that  the  above  definition  docs 

not  ensure  that  the  geodesies  p(6)   within  cr  for  9  €  (61,62)   are  actually  shortest  paths. 

However,  the  algorithm  to  be  described  below  will  make  sure  that  each  slice  a  that  it  will 

construct  will  correspond  to  some  peel  p  in  the  manner  just  described,  so  that  at  least  one 

geodesic  path  passing  through  the  sequence  of  edges  defined  by  ct  is  indeed  a  shortest  path. 

The  slice  tree  that  our  algorithm  will  construct  is  defined  as  follows.  Each  node  of  T  except 
for  its  root  is  a  slice.  The  root  is  a  dummy  slice  and  its  sons  are  the  slices  cr(6/,6/+ij/"o), 
i  =  l,  .  .  .  ,s,  where  /q  is  the  initial  face  of  K  containing  X  and  where  9i,  ...  ,8,  are  the 
orientations  of  the  segments  connecting  X  to  the  vertices  of /q,  ordered  in  angular  order  about  X. 
Let  a  =  a'(ei,02/)  be  a  node  of  T.  Its  sons  arc  obtained  by  extcixJing  cr  one  face  past  /. 
Specifically,  let  <  be  the  terminal  edge  of  cr,  i.e.  the  edge  of/  containing  all  the  terminal  points 
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Zq  for  6  d  (61,62),  and  let  Z'  be  the  face  of  K  adjacent  to/  at  tf.  If  the  geodesic  paths  in  a 
have  already  passed  through  f  (before  reaching/)  then  a  is  called  a  terminal  slice,  and  remains 
a  leaf  of  T.  Otherwise  we  extend  aU  paths  p(Q)  for  0  €  (61,62)  "^to/*,  and  let  T  denote  the 
set  of  points  at  which  these  paths  leave  /* .  Plainly  ""  ',i  a  connected  portion  of  the  boundary  of  /* 
and  thus  is  a  union  of  portions  Tj,  .  .  .  ,t„  of  edges  of/*  (actually  each  Ty  for  \<j<q  is  a  full 
edge).  For  each  l^j<q  let  »]iy  ^  (^u^z)  denote  the  orientation  of  the  geodesic  path  in  ct  which 
reaches  the  vertex  at  which  T;  and  Ty+j  meet.  Put  also  (JIq  =  6^  and  »{i,  =  62-  Then  slices  that 
are  candidates  for  the  sons  of  cr  in  T  are  the  slices  ct(»|i/_i,»|i,j/'),  i  =  l,  .  .  .  ,q  (see  Fig.  5.2). 


Bg.  5.2.  Extending  a  slice  p>ast  an  edge. 
The  problem  with  extending  the  tree  in  this  somewhat  uncontrolled  maimer  is  that  this  might 
cause  us  to  add  to  the  tree  slices  for  which  the  sequences  of  edges  of  K  which  they  define  may  be 
such  that  no  shortest  path  passes  through  it.  Note  that  the  number  of  true  slices  (i.e.  slices  for 
which  there  exists  at  least  one  shortest  path  passing  through  the  sequence  of  edges  which  they 
define)  is  0(/i  ),  in  virtue  of  Proposition  4.8,  but  that  we  have  no  obvious  way  to  estimate  the 
total  number  of  possible  "false"  slices  that  might  be  introduced  by  the  construction  in  the 
preceding  paragraph. 
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We  therefore  design  our  algorithm  so  that  it  ahvays  add  only  true  slices  to  the  tree  T.  More 
specifically,  the  algorithm  will  maintain  the  following  two  invariants: 

(a)  The  tree  constructed  by  the  algorithm  contains  only  true  slices  (in  the  above  sense). 

(b)  With  each  slice  a  =  o"  (6 1,6  2^0  ^^  associate  a  terminal  portion  ay  which  is  a  polygonal 
subregion  of  arClf.  The  terminal  portions  of  slices  in  the  tree  will  be  pairwi&e  disjoint  and  each 
terminal  slice  portion  will  contain  the  set  of  all  points  Z  in  /  for  w^ch  the  shortest  path  to  Z 
within  a  is  shorter  than  than  any  path  to  Z  contained  within  one  of  the  slices  present  in  T. 

To  maintain  the  invariant  (b),  the  algorithm  will  have  to  store  in  the  tree  that  it  constructs 
additional  information  concerning  the  terminal  portion  of  each  slice.  Terminal  slice  portions  that 
are  propter  subsets  of  their  corresponding  slice  portioiB  ad/  can  arise  due  to  the  overlapping  of 
such  a  afl/  with  other  slices  already  present  in  the  tree  (and  reaching  the  same  terminal  face). 
This  implies  that  as  the  tree  is  being  constructed,  these  trimmed  portions  may  be  trimmed  again 
and  again  as  new  slices  are  added  to  the  tree  until  they  reach  their  final  value  at  the  end  of  the 
algorithm.  We  will  show  below  that  proper  assembly  of  the  terminal  portions  of  the  slices  in  the 
final  tree  gives  us  all  the  peels  that  we  seek. 

Each  iteration  step  of  the  algorithm  adds  a  new  slice  to  the  tree.  Let  T/^  denote  the  slice- 
tree  after  the  t-th  iteration  of  the  algorithm.  For  each  starting  orientation  6  let  Z©  denote  the 
farthest  point  on  p(Q)  such  that  p(9^q)  is  whoDy  contained  within  a  slice  of  Tj^.  Define  the 
front  F(Ti^)  of  Ti^  to  be  the  set  of  all  points  Zg  which  lie  on  edges  of  ^.  Let  W{Ti^)  denote  a 
point  Z9  €  F(Ti^)  for  which  the  length  of  p(0^o)  is  smallest,  and  let  m(Tt)  denote  this 
smallest  length. 

Initially,  as  defined  above,  Tq  consists  of  all  the  slices  contained  in  the  face/o  containing  X, 
all  being  children  of  a  common  dummy  root. 

At  the  (t+l)-st  step,  the  algorithm  picks  Z^  =  W(Ti^),  and  obtains  the  slice  a  in  Tj^ 
containing  Z9.  Suppose  for  the  moment  that  Zq  lies  in  just  one  such  slice.  Let  /  be  the  terminal 
face  of  a,  let  c  be  the  edge  of  K  containing  Z^,  and  let/  be  the  other  face  of  AT  adjacent  to/  at 
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Let  t/9  denote  the  farthest  point  along  p (6)0/'  (lying  on  another  edge  of/*),  and  let  ct'  be  the 
slice  which  extends  a  past  e  and  which  contains  p{Q,Uq)  (again  we  suppose  for  the  moment  that 
only  one  such  slice  exists).  Then  we  have  '     -  . 

Lemma  5.1:  a'  is  a  true  slice,  in  the  sense  that  there  exists  at  least  one  shortest  path  passing 

through  the  sequence  of  edges  whidi  ct'  defines. 

...      ., ..  ^ 

Proof:  What  we  have  to  show  is  the  existence  of  a  starting  orientation  9'  within  the  range  of 
starting  orientations  of  ct'  and  a  point  U  ^  p{Q')(^f  for  whichp(6',l/)  is  a  shortest  path  to 
U.  We  take  9'  to  be  9,  and  f/  to  be  a  point  in  p(9)r\f  sufficiently  near  Zq.  Suppose  to  the 
contrary  that  p{Q,U)  is  not  one  of  the  shortest  paths  to  U.  Let  i|»  be  the  starting  orientation  of 
the  shortest  path  to  U.  Two  cases  can  arise: 

(i)  The  point  U  belongs  to  some  slice  in  Tj^.  But  if  this  is  the  case  for  all  points  U  sufficiently 
dose  to  Zq,  then  Zq  must  be  a  ridge  point  and  lie  on  the  boundary  of  more  than  one  slice  in  Tj^, 
contrary  to  assumption. 

(ii)  Hence  U  does  not  belong  to  any  slice  in  Tj^,  so  that  p(}\i)  must  reach  the  front  FiTj^)  at 
some  point  V  before  reaching  U.  But  then  the  length  of  p(ij/,V)  is  at  least  that  of  p(9^),  so 
that  if  U  is  chosen  sufficiently  near  Z^  the  path  p(»i<,t/)  cannot  be  a  shortest  path  to  U.  (To 
show  this  note  that  the  additional  length  of  p(9)  between  Zq  and  U  will  be  strictly  smaller  than 
the  length  of  p(«j')  between  V  and  U  if  one  chooses  U  sufficiently  near  Zg,  for  otherwise  an 
impossible  situation  as  in  case  (i)  above  would  occur.)  This  proves  the  lemma.  Q.E.D. 

The  cases  in  which  Zq  belongs  to  two  slices  in  Tjt,  or  in  which  there  are  two  slices  which 
extend  ct  and  contain  p{Q,Uq)  deserves  a  slightly  modified  treatment.  It  can  be  shown  that, 
except  for  some  rare  cases  which  can  be  detected  by  a  pxirely  local  analysis  of  the  structure  of 
geodesic  paths  near  Z^,  one  can  extend  either  of  the  slices  containing  Z^  in  the  first  case,  or  use 
any  of  the  extended  slices  obtained  in  the  seamd  case,  and  Lemma  5.1  will  still  hold  for  this 
extended  new  slice,  although  its  proof  will  have  to  be  somtwhat  modified. 
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The  algorithm  v\ill  then  add  the  slice  a'  as  a  son  of  ct  to  the  slice  tree,  thus  maintaining 
property  (a).  To  maintain  (b)  the  algorithm  will  have  to  trim  the  terminal  portion  of  ct'  to  keep  it 
disjoint  from  the  terminal  portions  of  other  slices  reaching  the  same  face,  and  possibly  also  trim 
these  other  terminal  portions. 

To  understand  how  such  a  trimming  is  to  be  done,  suppose  first  that  we  are  given  just  two 
slices,  both  reaching  the  same  terminal  face  /  and  overlappring  one  another  on  /.  Then  the 
following  basic  procedure  is  applicable. 

Slice  Trimming  Procedure:  Let  ct^  =  tr(9i,0'i/),  02  =  <^(02»®'2x/)  ^  ^o  *^ce*  with  the 
same  terminal  face/.  Apply  the  planar  unfolding  procedure  to  cr,  tr'  to  obtain  planar  triangular 
layouts  of  these  two  slices.  Move  these  two  layouts  in  the  plane  so  that  the  two  copies  of  the  face 
/  in  these  layouts  coincide  with  the  standard  plane  representation  of  /.  Let  X^, ,  X^  be  the 
positions  of  the  point  X  in  these  two  new  layouts,  and  let  /  be  the  perpendicular  bisector  of 
^a{^(T2-  ^^/l2  (resp.  /21)  denote  the  portion  of/  lying  on  the  X<,  -side  (resp.  the  Xg.  -side)  of 
/.  We  then  replace  o-jn/  (resp.  or2n/)  by  CTifl/  —  aino-2n/2i  (resp- 
CT2n/  -  CTinCT2n/i2).  see  Fig.  5.3. 
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Fig.  5.3.  The  slice  trimming  procedure. 
In  other  words,  for  points  Z  which  can  be  reached  by  two  distinct  geodesic  patlis,  lying  in 

two  different  slices,  the  slice  trimming  procedure  detennincs  which  of  these  two  paths  is  shorter, 

and  removes  the  point  Z  from  the  other  slice. 

Let  us  consider  in  more  detail  the  computational  aspects  of  this  procedure.  As  the  slice  tree 
is  being  built,  we  can  store  with  each  slice  a  the  position  X^  of  the  point  X  in  the  planar  layout 
of  a  in  which/  lies  in  its  standard  plane  position  (X^  will  be  represented  by  quantities  a^.,  Q„  as 
in  the  preceding  section,  which  can  be  easily  updated  as  we  extend  ct  to  an  adjfkcent  face). 

Given  two  slices  ctj  =  a(9i,0'ii/)  and  0*2  =  0'(92»^'2»/)  niceting  at  the  same  terminal 
face  /,  identify  cti  with  its  specific  planar  unfolding  used  in  the  slice  trimming  procedure  above, 
j  =  1,2.  Then  it  is  easily  seen  (cf.  Fig.  5.3)  that  a^  —  aiC\cr2C\f2[  is  a  bounded  polygonal 
region  which  is  star  shaped  with  respect  to  X^j^,  and  which  is  bounded  by  at  most  six  segments, 
including  the  two  rays  emanating  from  Xfj^.  We  regard  the  boundary  of  this  region,  ciduding  the 
two  rays  from  Xg.  ,  as  the  graph  of  a  function  expressing  the  distance  p  from  Xg.  as  a  function  of 
the  orientation  9  ^  (61, 6 '1)  and  write  it  as  Pa^jjX^). 
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Now  suppose  that  at  a  certain  time  during  its  construction,  the  tree  T  contains  slices 
a],  ...  ,a,  all  reaching  the  same  terminal  iac&  f.  Then  to  enforce  property  (b),  the  terminal 
portion  of  each  a,  will  have  to  be  replaced  by  v(CT,)n/,  where 

(*)  vM  =  ^i  -  U^ir^^Pfji 

Note  that  v(a,)  is  a  star-shaped  region  (with  respect  to  X„^)  whose  boundary,  except  for  the  two 
rays  emanating  from  X^. ,  is  represented  by  the  function 

p,/e)  =  min  p,,^/e) ,        9  €  (e„e',) 

In  other  words,  p^  is  the  lower  envelope"  of  if— 1  polygonal  lines,  each  consisting  of  at  most 
four  segments,  as  these  lines  are  viewed  from  the  point  X<j^. 

We  claim  that  the  following  property  holds: 

Lemma  5.2:  J£  p^j^  consists  of  f,  segments,  then  2^,  =  0(t). 

i 

Proof:  To  prove  this  claim,  we  proceed  in  a  manner  quite  similar  to  the  estimate  of  the  size  of 
Voronoi  diagrams  used  in  [Sti].  That  is,  we  define  an  undirected  graph  G  whose  nodes  are  the 
points  Xg.^,  and  in  which  an  edge  connects  X^  to  X^.  if  the  trimmed  portions  of  a,  and  (jj  have 
an  edge  in  common,  and  if  this  edge  is  not  contained  in  any  of  the  four  rays  emanating  from  X^.^ 
and  Xfj  .  It  is  easy  to  show  that  G  is  a  planar  graph  by  mapping  each  edge  (X^.^^- )  of  G  to  the 
union  of  two  segments  connecting  X^.  and  X„  to  a  point  on  the  corresponding  common  edge. 
Moreover,  the  graph  G  has  no  multiple  edges,  a  fact  which  can  be  established  as  in  the  proof  of 
Lemma  4.5.  It  therefore  follows  by  Eular's  formula  that  G  has  0{t)  edges.  Fmally,  each  of  the 
edges  obtained  by  the  trimming  process  which  is  not  recorded  in  G  lies  on  a  ray  bounding  some 
slice  CTj,  and  plainly  there  can  be  at  most  0{t)  such  edges.  These  observations  establish  our  daim. 
Q.E.D. 

Now  by  the  results  of  the  preceding  section,  /  :S  n,  since  the  total  number  of  "tnie"  slices 
reaching  the  same  face  is  at  most  n.  This  means  that  the  total  number  of  edges  separating 
trimmed  terminal  portions  of  slices  within  a  given  face  is  at  most  0{n)  during  each  stage  of  the 
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algorithm.  Moreover,  as  we  add  a  new  slice  a  to  the  tree  we  need  to  update  the  functions  p^.^  for 
all  slices  ct,-  reaching  the  same  terminal  face  /  as  O".  Suppose  as  before  that  Pg.  consists  of  f, 
segments.  Then  its  updated  value  due  to  the  appearance  of  a  can  be  computed  in  a 
straightforward  manner  using  (*)  in  time  0(f,),  so  that  updating  of  all  these  functions  can  be 
accomplished  in  time  0{n).  Moreover,  each  new  edge  added  to  any  of  these  functions,  and  only 
such  edges,  will  also  appear  in  the  graph  of  the  new  function  Pg.,  and  it  is  an  easy  matter  to 
assemble  all  these  edges  and  thus  obtain  p^  in  time  0(n). 

Remark:  Lemma  4.6  implies  that  the  final  trimmed  value  of  a  slice  is  convex  (when  properly 
unfolded).  The  procedure  sketched  above  might  however  temporarily  lead  to  nonconvex  trimmed 
slices.  We  do  not  know  whether  nonconvex  intermediate  slices  can  actually  be  obtained. 

We  therefore  conclude  that  one  can  maintain  property  (b)  in  time  0(n)  per  each  step  of  the 
algorithm.  '    ■ 

The  algorithm  also  needs  to  update  the  values  of  W{Ti^)  and  m{Tii)  in  view  of  the  addition 
of  ct'  to  the  tree.  Note  that  the  new  front  of  the  tree  is  obtained  from  the  previous  front  by 
deletion  of  the  "entering"  segment  of  cr',  by  addition  of  the  terminal  edge  of  ct',  and  by  possibly 
shortening  other  edges  of  the  front  due  to  the  trimming  procedure  described  above.  Note  that  the 
deletion  of  the  entering  edge  of  ct  '  will  in  general  have  split  (a  portion  of)  the  terminal  edge  of  ct 
into  two  subsegraents  which  still  belong  to  the  front,  and  that  similar  splits  or  replacements  of 
edges  will  result  from  the  trinuning  process.  This  suggests  that  we  represent  the  front  of  the  tree 
as  a  union  of  subsegments  of  edges  of  K.  For  each  segment  e  in  the  front  we  can  easily  compute 
the  point  W{e)  on  e  for  which  the  geodesic  path  to  W{e)  through  the  slice  bounded  by  «  is 
shorter  than  that  to  any  other  point  on  e,  and  also  the  length  fn{e)  of  this  path.  We  then 
maintain  a  priority  queue  containing  all  the  relevant  segments  e  constituting  the  front,  ordered  by 
the  value  of  m{e).  At  each  addition  of  a  new  slice  ct'  to  the  tree  we  delete  from  the  priority 
queue  the  terminal  segment  of  the  ancestor  ct  of  ct',  add  back  the  two  subsegments  of  this 
segment  still  in  the  front,  and  add  the  terminal  edge  of  ct'  .  Since  the  trimming  step  can  result  in  a 
new  terminal  slice  portion  having  no  terminal  edge,  and  since  the  terminal  edges  of  existing  slices 
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may  also  be  triramed  by  that  procedure,  the  priority  queue  of  front  edges  will  have  to  be  updated 
in  an  appropriate  manner.  Since  the  incremental  trimming  procedure  described  above  will  produce 
at  each  step  only  0{n)  slice-bounding  segments,  it  follows  that  at  most  0(rt)  updates  of  the 
priority  queue  will  be  required  at  each  stage.  Once  these  updates  are  performed,  the  updated 
values  W{Ti()  and  m{Tj^)  are  easily  available  and  can  be  retrieved  in  the  next  iteration  of  the 
algorithm. 

The  algorithm  terminates  when  the  priority  queue  representing  the  front  of  the  tree  becomes 
empty,  i.e.  when  the  front  itself  becomes  emptj'.  ■  t;  ■    .-.  .li  t       ■<-"  ; 

When  this  happens  we  still  need  a  final  phase  that  will  construct  the  peels  of  K  from  the 
slice  tree.  To  this  end  let  a  =  O'(0i,62»iO  ^  ^  ^^  o^  ^-  ^^  ^"^  trimmed  value  of  ct  is 
defined  as  .:.-,  .,       j.  ^  .  jisr. :..'.  •'  .  . 

ct' 

where  ct'  ranges  over  all  slices  on  the  path  in  J  to  ct,  and  where  t(ct')  denotes  the  fmal  trimmed 

terminal  portion  of  ct'.  That  is,  we  collect  all  trimmed  terminal  portions  of  slices  lying  along  the 
path  to  CT  in  r,  but  restrict  each  such  portion  to  the  wedge  between  the  two  starting  orientations 
defining  ct. 

Note  that  \i((t)  need  not  be  a  full  peel.  In  fact,  a  necessary  and  sufficient  condition  for 
fJL(CT)  to  be  a  peel  is  that  the  two  bounding  geodesies  p(9i)  andp(Q2)  are  not  trimmed,  and  still 
reach  vertices  of  K  within  (the  closure  of)  jjl(ct).  K,  say,  p(Q\)  has  been  trimmed,  let  ct^  be  the 
slice  whose  range  of  starting  orientations  is  adjacent  to  that  of  ct  at  Sj.  Tben  it  is  easy  to  see  that 
the  peel  containing  ft(CT)  also  contains  }Jl(cti).  These  observations  makes  it  dear  that  we  can 
construct  all  peels  by  a  depth-first  traversal  of  T,  visiting  sons  of  a  slice  in,  say,  counterclockwise 
order  of  their  ranges  of  starting  orientations. 

To  estimate  the  time  required  by  the  algorithm  we  note  that  the  maintenance  of  property  (b) 
is  the  costliest  part  of  the  algorithm,  in  which  the  updating  of  the  priority  queue  representing  the 
front  of  T  may  require  0(n  log  n)  time  for  each  step  (the  trimming  procedure  itself  requiring 


-39- 

only  0(n)  time).    Snce  the  algorithm  adds  at  most  0(n^)  slices  to  the  tree,  it  follows  that  the 
algorithm  will  run  in  time  0{n    log  n). 

The  correctness  of  the  algorithm  follows  from  the  following  considerations.  First  note  that 
the  peels  as  constructed  by  the  algorithm  are  pairwise  disjoint  by  construction.  We  also  claim  that 
they  cover  the  whole  surface  o(  K.  To  siiow  this,  it  suffices  to  prove  that  every  true  slice  is 
constnicted  by  the  algorithm.  Indeed,  if  this  latter  property  is  known  to  hold,  then  for  each  point 
Z  on  the  surface  of  K  let  pC^v-^o)  '^'"^  '*-  s-iOilwt  path  to  Z.  The  sequence  of  edges  and  faces 
through  which  this  path  passes  defines  a  ti-ue  slice  a  which  will  appear  in  the  final  tree  T.  It  is 
easy  to  verify  that  Z  wiU  belong  to  the  final  trimmed  terminal  portion  of  cr. 

To  see  that  all  true  slices  are  coi'structe'J  by  the  algorithm,  we  note  that  once  the  algorithm  has 
added  a  slice  to  the  tree  with  a  nonempty  terminal  edge  which  cxmoects  that  slice  to  another  true 
slice,  then  this  edge  will  not  be  wholly  deleted  by  the  trimming  procedure,  and  eventually  the 
following  slice  will  also  be  picked  up  by  the  algorithm  and  added  to  the  tree.  This  implies  in  a 
straightforward  inductive  manner  that  all  true  slices  are  added  to  the  tree. 

Hence  the  peels  produced  by  the  algorithm  are  pairwise  disjoint  and  cover  the  whole  surface 
of  K.  It  is  now  a  simple  matter  to  prove,  arguing  as  in  the  proof  of  Lemma  5.1,  that  each  such 
peel  is  one  of  the  peels  defined  in  the  preceding  section,  and  thus  to  conclude  the  proof  of 
correctness  of  the  algorithm. 

In  view  of  the  discussion  at  the  preceding  section,  we  thus  have  the  following  summary  theorem. 

Theorem  5.2:  Given  a  convex  polyhedron  K  with  n  vertices  and  a  point  X  on  its  surface,  one  can 
preprocess  AT  by  a  procedure  which  runs  in  0{n  log  n)  time.  This  procedure  produces  a  data 
structure  of  size  0(n^,  with  the  aid  of  which  one  can  find  in  0{n)  time  the  shortest  path  along 
the  surface  of  K  from  X  to  any  other  specified  point. 

Remark:  It  seems  quite  likely  that  the  algorithm  developed  in  this  section  is  not  optimal,  as  it 
requires  0(n  log  n)  time  to  construct  a  quadratic  data  structure  (and  then  search  the  structure  in 
only  linear  time  to  find  a  shortest  path).  After  the  original  submission  of  this  papa-,  Mount  [Mo] 
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has  recently  proposed  an  improved  approach,  in  which  the  data  structure  only  maintains  points  of 
intersections  of  slices  with  the  edges  of  K  (rather  than  with  its  faces,  as  done  here).  Thus  his  data 
structure  is  a  collection  of  disjoint  intervals  on  the  edges  of  K,  rather  than  a  collection  of  dbjoint 
polygons  on  the  faces  of  K,  making  it  much  easier  to  maintain  this  structure,  and  thereby  reducing 
the  running  time  to  C>(n  log  n). 
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